“Calhoun 


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 


DSpace Repository 


Theses and Dissertations 


1973 


1. Thesis and Dissertation Collection, all items 


Hydraulic ram studies. 


Soper, Wesley Richard. 


Monterey, California. Naval Postgraduate School 


http://ndl.handle.net/10945/16644 


Downloaded from NPS Archive: Calhoun 


agitth | DUDLEY 


NN \ KNOX 


LIBRARY 
http://www.nps.edu/library 






Calhoun is the Naval Postgraduate School's public access digital repository for 

research materials and institutional publications created by the NPS community. 

Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 
appointed — and published — scholarly author. 


Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





HYDRAULIC RAM STUDIES, 


— 

ae. «bias 
/ VS, 
| - « 





NAVAL POSTGRADUATE SCHOOL 


Monterey, California 





THESIS 





HYDRAULIC RAM STUDIES 
by 


Wesley Richard Soper 


Thesis Advisor: Hee ige “owe r 


December 1973 


Approved for public release; distribution unluncted. 





Hydraulic Ram Studies 


by 
Wesley Richard,Soper 


Lieutenant, United States Navy 
B.S., United States Naval Academy, 1967 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN AERONAUTICAL ENGINEERING 


from the 


NAVAL POSTGRADUATE SCHOOL 
December 1973 








Library 
Naval Posters Ccnoot 
Monterey, Califo 33940 


ABSTRACT 


Hydraulic ram concerns catastrophic failure of liquid 
filled tanks impacted by high speed projectiles. It is 
divided into four phases: entry/shock, pressure field, 
cavity, and free surface. The entry wall phase was investi- 
gated in detail. 

A formula to predict energy losses of a bullet penetrating 
an aluminum plate is presented, and experimental data were 
obatined and compared. 

A fiber optic transducer was designed to measure generated 


pressures in the shock and cavity phases. 
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LI. INTRODUCTION 


Since the days aircraft first began to fly in a hostile 
environment, the achievement of non-nuclear aircraft sur- 
vivability has been the subject of extensive research. The 
ability of today's expensive air-to-air and air-to-ground 
weapons delivery systems to withstand hits by small arms 
fire and by warhead fragments is of paramount importance in 
this era of decreasing defense budgets. 

An aircraft component which is easily hit and most 
susceptible to damage is the fuel tank. A damaged fuel cell 
can lead to aircraft loss in several ways, some of which are 
fire, explosion, fuel starvation, or fuel cell failure. 

Catastrophic failure of a fuel tank due to impact by a 
high speed projectile is known as hydraulic ram. This is 
actually a misnomer, since the failure mechanism consists of 
four distinct phases rather than one quasi-static compression 
type of state. [Ref. 1] 

The first phase, the entry/shock phase, is caused by the 
projectile compressing the fluid near the entry point to 
extremely high pressures. These pressures can range in the 
hundreds of thousands of pounds per square inch and cause a 
shock wave to move radially outward from the jorelatoe nei. Cichmie ae 
However, the shock attenuates rapidly and, in general, is 
too weak to cause damage to the side and rear walls of the 


fuel tank; but it can cause outward bending of the wall in 





the vicinity of the entry point and start cracks which can 
lead to eventual destruction of this wall. 

At any given point in the tank a (relatively) slow 
pressure pulse begins when the attenuated shock wave reaches 
that point and continues with the passage of the projectile 
through the liquid. This phase is called the pressure field 
phase. Its pressures range in the hundreds of pounds per 
square inch and last for several milliseconds. 

The third phase of hydraulic ram is called the cavity 
phase. As the projectile eee through the fluid it must 
displace the fluid from its path. A large radial velocity 
is thus imparted to the fluid, and, as a result, the fluid 
moves away from the projectile path creating a cavity. 
Existing fluid pressures (gravity and pressures caused by 
tank boundaries) eventually halt cavity Prcuctoana cause the 
onset of collapse. Inside the cavity are fluid vapors and 
trapped air which entered during wall penetration. As the 
cavity collapses this air is compressed to a high pressure 
resulting in re-expansion of the cavity. Expansion and 
collapse repeat several times before all the stored energy 
is dissipated. Pressures during this phase are of the order 
of one thousand pounds per square inch. 

These three phases combine to impart a large velocity 
to the upper free surface of the liquid, if ullage is present, 
which in turn imparts an impulsive load on the top of the 
tank. This fourth phase of hydraulic ram is known as the 


free-surface phase. 





The structural response to hydraulic ram is brought 
about by the transfer of energy (projectile to fluid, fluid 
to walls, etc.). The failure mode is a function of how this 
energy is distributed among the four phases. In this study 
the projectile energy loss during entry wall penetration was 
investigated in detail and compared with experimental 
results. There are two areas of bullet penetration studies: 
thosé assuming a plate of finite thickness which the projec- 
tile passes through and those concerning a plate of. thick- 
ness greater than the penetration depth. The effect of 
projectile shape on penetration distance and energy loss by 
a projectile passing through a finite thickness plate were 
Ronainened in this study. In order to be able to predict 
energy transferred to the fluid by the projectile during 
the shock phase and later during the cavity phase, a know- 
ledge of energy losses during entry wall penetration is 
Sssicibial « 

Knowledge of pressures during the three final phases of 
hydraulic ram is important because the major damage occurs — 
auring =Gnem. sstiber Optic pressure’ transiducers were designed 
in order to experimentally measure these pressures. The 
magnitude and duration of pressures during the shock and 
cavity phases were examined, and transducers were designed 


which are capable of accurately measuring these characteristics. 





IIT. ANALYTICAL PENETRATION STUDIES 


In order to predict the phenomena taking place as a 
projectile impacts and passes through a plate, an analytical 
study of the entry wall phase of hydraulic ram is presented. 

An essential assumption of this study of projectile 
impact, and the assumption made by Dunn in Ref. 2, is that 
"Chere exists in each unconstrained material body at constant, 
uniform initial temperature a dynamic yield strength, Ov. 
that is a constant of the material." Therefore there is one 
well-defined material constant which can be used to predict 
penetration data. 

This dynamic yield strength, when multiplied by the area 
of the projectile, A(x), which resists the motion, gives an 
estimate of the force retarding the bullet. An additional 
assumption is that the dynamic yield strength of the 
projectile is much greater than that of the target so that 
no substantial deformation of the projectile occurs. For 
these assumptions the retarding force is given by 

F = =e A(x) es) 
where oy is the yield. strength of the target. Thus the 


Ah 
projectile equation of motion is 


0 
re Yen 
Xx + = ACY) = Q C23) 


where x is the acceleration of the projectile mass m. One 





Salient feature of equation (2) is that the force retarding 
the projectile is only a function of x. It is assumed that 
there is no plastic deformation of the target material and 
no significant frictional force during penetration. 

Equation (2) was used for two types of penetration solu- 
tions: those concerning impact into a plate of thickness 
greater than penetration depth, and those concerning the 
passage of projectiles through thin plates. 

For the first case, multiplying (2) by x and integrating 
yielded 


P 
(=2M)* (x° — v,°) = - f a(x) ax (3) 


= 
where G is the shear modulus of the plate and approximately 
equal to eres (Ref. 2], Vy is the impact velocity, and P 
the penetration depth. 

Equation (3) was solved numerically by a computer 
program. The program, for a given impact velocity, deter- 
mined the bullet trajectory parameter V as a function of 
incremental values of x until zero velocity was reached at 
the penetration depth. Figure 1 shows penetration depth vs. 
impact velocity for forty-five grain projectiles of three 
shapes. As expected, projectiles with a smaller cross- 
sectional area gradient penetrated further for a given impact 
Verocicy. 

The second penetration problem consicered was that for a 


given plate thickness. Unlike the first case, where the 


10 





effective retarding force area distribution A(x) increased 
until it reached a constant maximum value, A(x) for a thin 
plate of given thickness increased only until the nose of the 
projectile penetrated the backside of the plate. At that 
time A(x) decreased due to projectile protrusion from the 
backside of the plate. Figure 2 shows two stages of pro- 
jJectile penetration. It should be noted that the x- 
coordinate moves with the projectile so that x is defined 

as the distance from the nose of the bullet to the front 
side of the plate. In Figure 2a the bullet has partially 
penetrated the plate. The effective retarding area is shown 
In Figure 2b the nose of the projectile has 


as Arr: 
penetrated the back side of the plate and A. is the 


ff 

annular area shown. Thus, as the maximum cross-sectional 
area reaches the backside of the plate, the area resisting 
motion goes to zero. It is at this point that penetration 
is complete. The minimum impact velocity which a projectile 
must possess to ensure its reaching this point is defined as 
Poe valwisvic Limit. 

The work done on the projectile during its deceleration 
in the plate is equal to the change in kinetic energy during 


penetration. An integration of the total retarding force 


during penetration yeilds 


W=oa A t (4) 


where tvis' the thickness of the plate. It is seen that the 


volume swept out by the projectile determines the penetration 


an 





energy loss, and the energy loss is totally independent of 
projectile shape. Alternately stated, this result predicts 


that, for a projectile of given mass and A. the ballistic 


ax? 
limit is independent of projectile shape. 
With this in mind, equation (2) was solved assuming a 
constant cross-sectional area, that is, assuming a cylindri- 
cally shaped projectile. The equation of motion was inte- 
grated twice with appropriate initial conditions substituted 


to yield: 


ea te TG (5) 


where T is time and A is the projectile cross-sectional area. 
For a projectile with a velocity of V; impacting a plate 


of thickness t and exiting with a velocity of Ve> from (5) 


above, 
2 
4 GA TS 
oS ieee 15m 2 (6) 
and 
Ve i ee (7) 


folvinems(;)) for T and substituting into (6) and®solving for 


‘- yields 


wo~ pee Bebe F 


We 
e alt 15m (8) 


memeene Dalilieuuwc, limit. the exae velocity is zero, hence 


equation (8) was solved for V5 (denoted as V, to designate 


ie 





ballistic limit): 


ae CA 
Vo = 15m t J 


We (9) 
A relationship between the impact velocity, the exit 
velocity, and the ballistic limit can be found by substituting 
(9) into (8) and solving for V,: | 
V 


= wig 
Ve es [ 7 -1] 


LC (10) 


Equation (10) shows, for impact velocities much higher 
than the ballistic limit, that the exit velocity approaches 
the impact velocity. 

Equation (10) may also o used to find an estimate of 
the percentage of energy loss for a projectile which pene- 
trates a finite thickness plate: 

Z GA 


Se * a 
Lone [ E, |] 100 = [ ism ¢ eule a |] LOGy (11) 


a 


%E 


where Eo is the projectile kinetic energy when traveling at 


Ebemballistic JNimityeandPE, isithe energy eat impact velocity- 
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a. Early Stage 


A A(x) A(x-t) 





Arr = A(x) -— A(x-t) 
b. Protrusion Stage 


FIGURE 2. Stages of Projectile Penetration 
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IIL. EXPERIMENTAL STUDY OF ENTRY WALL PENETRATION 


Figure 3 shows the basic components of the ballistic 
range used in the experiments. They were the rifle and its 
mount, chronograph screens and counters, aluminum test plate 
and holder, shadowgraph station, and bullet catcher. 

The rifle mounting system, shown in Figure 4, consisted 
of the rifle mount and an adjustable rigid rifle aeune 
stand. The system was designed to hold rifles of .22, .222, 
and .30 caliber. In this are, a .222 Remington was used. 
The bullets utilized were of a mass of forty-five grains and 
of two shapes, hornet and spitzer. Projectile impact 
velocities ranged from 2400 to 3000 feet per second. 

The velocity was measured by three chronograph screens, 
AVTRON No. A914T333, mounted six feet SOeRe - Each screen 
had a five-volt D.C. signal across it. When the screen was 
broken by a bullet, a pulse was sent to a counter. The first 
sensor started the first counter. The second sensor simul- 
taneously stopped the first counter and started the second. 
The third sensor stopped the second counter. Thus the 
average velocity between counters was obtained, allowing a 
prediction of the impact velocity at the plate. A typical 
chronograph station is shown in Figure 5. Also shown in 
this figure is a shadowgraph apparatus which was not utilized 
dUmgumemenis phase of testing.» The counters wsed were 
Monsanto 101B 1MHz. The trigger error was less than © 0.3% 


of one period for sine waves with signal-to-noise ratio of 


40 db or better. 
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The empty fuel tank entry wall was simulated by a three- 
inch by four-inch 7075-T6 aluminum plate. Three thicknesses, 
0.05 inches, 0.09 inches, and 0.16 inches were tested. 

Karly in the test phase it was determined experimentally 
that the size and method of supporting the test plate had no 
significant effect on the exit velocity of the projectile. 

In addition to stopping the 4 Bewad counter, sensor three 
was used to trigger the delay unit which generated a pulse 
which could be time delayed from 190 useconds to 1700 
useconds. The delay unit circuitry is shown in Figure 6. 

The delayed pulse was then used to trigger a spark source 
just after bullet penetration of the plate. The light from 
the spark source was directed to a film plate by a concave 
mirror on the opposite wall. The mirror collimated the light 
which then passed behind the test plate and onto the film. 

In this manner, shadowgraphs were taken to examine the shape 
of the bullet after it had passed through the plate. The 
shadowgraphs also showed the resulting flow field behind the 
plate. Figure 7 shows the spark source and mirror as they 

WC TeasS C lea Dr. 

The chronograph screens mounted directly behind the test 
plate were used to measure the velocity of the bullet as it 
exited the plate. Figure 8 is a photograph of this apparatus. 

The bullet’ caltcher is shown in Figure®9. The sides, top,-> 
and bottom were 1/2-inch aluminum. The three-foot square 
front had three 3/8-inch plywood baffles inserted into it. 


The back of the catcher was made of two plates of 3/8-inch 


IP ge 





steel mounted at an angle of 45 degrees to the path of the 
bullet. This ensured bullet deflection downward into a 


layer of sand, reducing the likelihood of ricochets. 
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Ballistic Range Rifle Mounting System 


FIGURE 4. 
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IV. FIBER OPTIC PRESSURE TRANSDUCER DESIGN 


A. INTRODUCTION 

An accurate means of measuring pressures generated 
during all phases of hydraulic ram is needed to assess the 
validity of any analytical studies. The most popular types 
of pressure measuring SOHLCRE in use today are the semi- 
conductor strain gauge and the piezoelectric transducer. 
Both costly and fragile, these devices are also sensitive 
to temperature changes and contamination from small dust 
particles. A fiber optic pressure transducer was designed 
in order to overcome these problems. Two optical glass 
fiber bundles were used, one to transmit light to the inner 
surface of the transducer diaphragm and one to collect and 
transmit reflected light to a photomultiplier tube. This 
signal output was then amplified and could be measured from 
a photograph of the resulting oscilloscope trace. The 
intensity of the output light, being a function of displace- 
ment between the fiber bundle and the diaphragm, was thus 
made a measure of pressure on the diaphragm. A possible 
arrangement of the pressure measuring apparatus is shown in 
Figure 10. 

The fiber optic transducer is in reality a proximity 
gauge. Because pressure on the diaphragm causes it to 
deform slightly and decrease the distance between the 
Glaphragm surface and the fiber ends, output can be read in 


terms of pressure increase rather than distance decrease. 
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Figure lla, obtained from Ref. 6, shows two fibers, one 
emitting light and one collecting light. As the displacement 
between the fiber ends and the diaphragm increases, the cone 
of light from the transmitting element illuminates an 
increasingly larger area of the diaphragm. The light 
reflected from the diaphragm becomes a cone of light, 
increasingly intense as the diaphragm moves. away from the 
fiber ends. Figure llb shows typical output of the 
receiving fibers vs. fiber-diaphragm displacement. The 
front side of the curve indicates linearity \Ibetween output 
and distance. As the distance increases, however, the cone 
of light from the transmitting fiber fills tthe entire 
diaphragm, and the curve reaches its peak. [Beyond this 
point, the intensity of light to the receiviung element 


decreases as the inverse of the distance sqmared [Ref. 6]. 


B. DESIGN CONSIDERATIONS 

Due to the great disparities between the peak pressures 
and the time scales of the pressure signals ‘tin the shock and 
cavity phases of hydraulic ram, two separate pressure trans- 
ducer designs are necessary. Figure 12 shows the basic 
geometry of both transducers. In addition to the optical 
fiber probe, the transducers consist of an outside tube and 
a diaphragm silver soldered or welded to the tube end. A 
sleeve is required inside the tubing to ensure orthogonal 


orlentation of the probe tip to the diaphragm. 
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Linear response of the pressure transducers i desired, 
and in order to achieve this the maximum deflection of the 
diaphragm is limited to less than half its thickness. This 
allows the light transmission to remain on the linear 
portion of the displacement curve of Figure ll. The resonant 
frequency of the diaphragm must be large enough to provide 
adequate frequency response to the transducer. Finally, the 
pressure transducer must be able to withstand the maximum 
pressures occurring during hydraulic ram and should be placed 
so that the Rnonee an is flush with the tank wall when 
measuring shock phase pressures. The transducer should be 
isolated from wall vibration, if sans therefore, it was 


insulated by a seal made of RTV. 


C. SHOCK PHASE DESIGN 

Figure 13 shows peak pressures (predicted by the strong 
shock theory of [4]) encountered at various distances from 
point of impact of a .30 caliber projectile penetrating the 
fluid at 3000 fps and an energy of 34,000 ft-lb. Measurable 
pressures are of duration of about 12 usec and a pressure 
transducer placed approximately 1.4 inches from the point of 
initial impact will experience a peak pressure of 7.8 x 10" 
psi, l2 psec after impact. Figure 13a indicates that the 
shock pressure signature can be approximated by the sawtooth 


fumetOnwshown in FPiguresl3b. The Fourier Expansion of the 


sawtooth is given by 
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n 
are < 
f(x) =7+ ) — cos kx + (1) sin kx } O<x<2m (12) 
kK=133,5.. tek 
k 
_ 2 _ (1) 
where ay om. and Dy ine a 


1 k 


For adequate frequency response of the transducer, it 
is desired that the amplitude of the high frequency signal 


content be at most ten percent of the fundamental amplitudes 


ay and Dy. For ay proportional to cae the smallest value of 


k that meets this requirement is 5. However, DL 
tional to k7! and thus the smallest value of k is ll. The 


is propor- 


high frequency content of the shock pressure signal is then 
approximately k/24 usec or 0.46 MHz. 

Requiring the transducer to be capable of measuring peak 
pressure of 100,000 psi, Figure 15, obtained from Ref. 3, 
shows that the minimum diaphragm thickness is 23 mils. 
According to Figure 16, all stainless steel diaphragms of 
sufficient diameter (greater than 0.160 inch to accommodate 
the .109-inch fiber bundle) have resonant frequencies near 
the operating frequency. A diaphragm 35 mils thick and 0.160 
inches in diameter has a resonant frequency of 0.5 MHz and 
hasea maximum defliection of just under 1 mil at 100,000 psi. 
Because the fiber optic system output resolution is approxi- 
mately 0.01 mils, only pressures above 1000 psi can be measured 


PO@eehies Dreposed,desiean. 
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D. CAVITY PHASE DESIGN 

The geometry of a transducer measuring pressures during 
the cavity phase must be different from that required to 
measure pressures during the shock phase due to lower pres- 
sures acting over longer periods of time during cavity 
collapse. 

Figure lta, obtained from Ref. 5, indicates, for a 0.30 
caliber projectile, maximum cavity collapse pressures were 
less than 1000 psi and each pulse lasted for approximately 
0.5 msec. The pressures rose and fell at nearly the same 
rate, allowing the pressure signature to be approximated by 
a triangular wave, as shown in Figure 1b. 

As in the previous analysis, the Fourier expansion of 


the triangular function ytelds: 


f(x)=5-5 } +5 cos kx 0<x<2n (13) 
TY) k=1,3,5 2K 
Again the amplitude of the high frequency components must be 
10% of the fundamental amplitude for adequate frequency 
response. The frequency response required by this pressure 
signal is then k/0.5 msec or .01 MHz. 

A stainless steel diaphragm 10 mils thick with a 
diameter of 0.190 inches has a resonant frequency of very 
near O7l MHz, and is capable of withstanding peak pressures 
of 3000 psi. It is suggested that this geometry be used 


for the cavity phase transducer. 
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V. RESULTS AND CONCLUSIONS 


One important assumption made in the analytical study was 
that the yield strength of the projectile was significantly 
higher than that of the target material. However, the pro- 
jJectile was copper jacketed, and the yield strength of copper 
is approximately 0.39 x 10° psi as compared with 0.25 x 10° 
psi for aluminum. Figure 17 shows clearly that, although 
the projectile and target yield strengths are of the same 
order of magnitude, projectile shape was not altered during 
plate penetration, except very close to the nose of the 
projectile. Thus it can be concluded that for hypervelocity 
impact phenomena, the assumption only need be made that the 
yield strengths of the projectile and target materials be 
of the same order of magnitude. 

iow diaiwercaim rcouly ulal iprOJ,CeCLI le esnapes 1S nov 
material in preterit penetration energy losses is also 
shown in the experimental results. Equation (9) was derived 
using that assumption. Figure 18 shows a plot of exit veloc- 
ity vs. impact velocity for 45-grain projectiles of two 
shapes penetrating aluminum plates of three thicknesses. 

The solid lines are plots of Equation (10), using Equation 
(9) as the definition for ballistic limit. Because both 
Shapes yielded similar data, no distinction between them was 
made. Data points agreed well with prediction leading to 
CiemeonclUsion tThateprojectile shape dods mot, in fact, 


significantly affect energy loss through a plate. 
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The theoretical energy loss from Equation (11) was 
plotted vs. impact velocity in Figure 19. Experimental 
data again agree closely aCe prediction, reinforcing the 
concept that projectile shape has no significant effects 
on energy loss through thin plates. 

It can be concluded, therefore, that a means is avail- 
able whereby energy losses through a wall of an unfilled 
fuel tank can be predicted accurately and simply. This is 
the first step in finding energy transferred during the 
hydraulic ram phenomenon. Effects of fluid foundation on 
the penetration characteristics of the walls are expected 
to be minimal, but further experimental study of this effect 
Should be undertaken. 

It should be noted that all experimentation and analyti- 
cal work in this study has been conducted using bodies of 
revolution that do not tumble. However, an application of 
the concept of the shape independence of energy loss can be 
made on Se case of fragment penetration. This representa- 
tion will only be applicable to early stages of hydraulic 
ram, aS projectile shape does have major effects on the size 
and intensity of the cavity phase. Aircraft damage and 
penetration characteristics of fragments (such as from SAM 
warheads) may be approximated by cylindrical bodies of equal 
mass and maximum cross-sectional area. It is expected that 
fragment energy loss during wall penetration will follow 


Equation (11). 
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VI. RECOMMENDATIONS FOR FURTHER STUDY 


Experimental studies in the area of hydraulic ram are 
continuing and further information will be gained. A few 
specific recommendations are presented as guidelines for 
future research. | 

The ballistic range is equipped to fire .30 caliber 
projectiles. This should be done in order to obtain data 
for projectiles of higher energies and larger damage radii. 
In addition, .30 caliber projectiles offer a greater 
variety of nose shapes which can add information to the 
previous data on bullet shape effects. 

Filled tanks should be penetrated by projectiles of 
various caliber and mass. This will yield information 
concerning the effect of fluid foundation on penetration 
of thin walls. 

Finally, it is recommended that the pressure transducers 
be installed and that experiments be conducted to confirm 
their effectiveness. If successful, this would give 
pressures as functions of time for the shock and cavity 
phases of hydraulic ram, a major step in understanding the 


phenomenon. 
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APPENDIX A 


Explanation of Computer Program 


This program computes the distance which a projectile of 
arbitrary shape and mass will penetrate into a semi-infinite 
metal plate. Inputs to this program are the shear modulus 
of the plate material, the mass of the bullet, its maximum 
radius, and its initial velocity. The cross-sectional area 
of the bullet at each .032 inch aft of the nose is also read 
in. Equation (3) is integrated trapezoidally using the 
bullet area distribution and velocity is tested until it 
becomes negative. The previous velocity is essentially zero, 
and the corresponding distance is output as the penetration 


depth. 
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